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An approach is proposed to analyze the azimuthal distribution of particles produced on the away 
side in heavy-ion collisions without background subtraction. Measures in terms of factorial moments 
are suggested that can suppress the statistical background, while giving clear distinction between 
one-jet and two-jet event structures on the away side. It is also possible to map the position and 
strength of the recoil jet to suitably chosen asymmetry moments. 
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An important area of investigation in relativistic 
heavy-ion collisions is the azimuthal angular distribution 
of hadrons produced on the opposite side of particles trig- 
gered at high pt 0, 0, 13 ■ That distribution reveals the 
properties of the medium effect on a hard recoil parton 
traversing the dense system formed by a central colli- 
sion. Several possibilities have been suggested for what 
may happen to the hard parton and what signatures may 
be registered on the away side 0, 0- ■ Our concern in 
this paper is not so much on the nature of the possi- 
ble signatures, but on how best to extract them from a 
sea of noisy background that is inherent in every nuclear 
collision. The conventional approach is to subtract the 
background. That is accomplished by first summing over 
all events subject to specific kinematical cuts. It is mean- 
ingless and impossible to make background subtraction 
event by event. However, summing over all events is a 
process in analysis that often degrades the signal, since 
the recoil parton may be absorbed by the medium, or 
may emerge on the other side as a reduced minijet, or 
may generate a shock wave. To enhance the signature 
what is needed is a measure that filters out the statisti- 
cal fluctuation in every event to the extent possible, and 
then let that measure be summed over all events. The 
aim of this paper is to propose such a measure and to 
demonstrate its effectiveness in the framework of a sim- 
ple model that simulates events with jets in the midst of 
a large statistical background. 

The measure that we propose is factorial moments. Be- 
fore going into the details, let us first give some historical 
background. The factorial moments were first used in the 
study of intermittency, which quantifies the self-similar 
behavior of branching processes in multiparticle produc- 
tion 0,0- It has also been suggested for critical behavior 
in heavy-ion collisions [jj. Although intermittency has 
been found in elementary collisions, its study in nuclear 
collisions has been plagued by a number of difficulties, 
among which is the problem of statistical fluctuations in 
large systems despite the theoretical virtues of the fac- 
torial moments. There are several differences between 
those intermittency studies and our problem here. First, 



we do not analyze fluctuations in rapidity, and do not 
consider a wide range of bin size to search for fractal be- 
havior. Second, we consider only a subset of rarer events 
selected by high pt triggers, and examine the away-side 
Acf> distribution that consists of far lower multiplicities 
of particles. Third, in intermittency studies the dynami- 
cal fluctuation is convoluted with statistical fluatuation, 
whereas our jet signal on the away side is additive rela- 
tive to the background. Having made these remarks to 
dissociate our work from whatever vestige there may be 
from the past, we proceed now to a description of our 
measure from the basics. 

The usual A(j> variable is defined relative to the trigger 
momentum in the transverse plane. For the convenience 
of our analysis on the away side, let us define ip — A0— tt. 
Consider an interval / around tp = 0, which, for definitc- 
ness, may be taken to be from —1.5 to +1.5, although 
the suitable range is an experimental decision. Let / be 
divided into M equal bins so that the bin size is 6 = I/M. 
In an event let rij denote the number of particles in jtYi 
bin. Define the factorial moment (FM) by 

1 M 

i=i 

and the normalized factorial moment (NFM) by 

p q = u/n, (2) 

for each event. It is important to recognize that Eq. 
(1) is used to determine the FM event by event, and 
the summation there may be regarded as the horizontal 
average. 

Let us now consider the hypothetical case where the 
fluctuation of rij from bin to bin can be described by 
a Poisson distribution as an ideal representation of the 
background 

n n 

P ri {n) = —e^, (3) 
n\ 

where n is the bin multiplicity that fluctuates and n is the 
average that depends on <5 and total event multiplicity N . 
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Then Eq. (JTJ can be written as 

E 



(n-?)l 
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If N is large compared to n, it can be approximated by in- 
finity in Eq. (@J; it then follows from Eq. © that f q — n q , 

so we have Fq Stat ^ = \ fo r any 5 in the ideal statisti- 
cal case. When N is finite but large, the corresponding 

f stat") 

Fq should be slightly less than 1. Furthermore, it has 
been pointed out to us that F q as defined in Eq. J2J i s 
independent of detector efficiency ^(| . 

In reality, F^ stat ^ [ s not exactly 1. We expect F q to 
fluctuate from event to event. Let us define in general 
the event-averaged NFM 
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where iVj is the multiplicity of the ith event and A/" ev t 
is the total number of events. The central question is 
whether {F q ) stays approximately at 1 for pure back- 
ground in the real data and becomes larger than I signif- 
icantly enough, when jets exist, so that (F q ) can be used 
as an effective measure to distinguish the signal from the 
noise. We stress that the (vertical) average over all events 
in Eq. (J5J is performed after the horizontal average is 
done in Eq. l(T]l. not the other way around. Thus if the 
answer to the above central question is in the affirma- 
tive, then one may regard F q as effectively suppressing 
the statistical background event-by-event. Since it is im- 
possible to make background subtraction event-by-event, 
to suppress the statistical contribution to a suitably cho- 
sen measure seems to be the best that one can do, and it 
should be done at the level of each event. 

To test the effectiveness of the measure suggested 
above, let us consider a simple model to simulate the 
tp distributions of events with jets in the presence of sig- 
nificant background. Let Ni particles be distributed ran- 
domly in the interval / with — f .5 < tp < 1.5. They will 
be regarded as the background of the ith event. For now, 
we limit ourselves to the simplest case of N% being a con- 
stant, but later generalize iV, to be Gaussian distributed. 
For a 1-jet event, we add a cluster of particles of multi- 
plicity, m = 5, bunched in a small interval in tp of extent 
e = 0.04. The cluster is randomly located in the interval 
— f < tp < +1. For a 2-jet event (to mimic a Mach cone 
structure) we put in two clusters of 5 particles each. We 
label the three cases by bg, lj and 2j, respectively, where 
bg denotes background only without jets. In Fig. 1 we 
show the distributions of (a) F2 and (b) F3 for the three 
cases, for Ni = 60, M = 30 and 7V cvt = 2000. Although 
the distributions have significant overlap, their average 
values (F q ) are distinctly separated, as shown in Fig. 2, 
for various bin numbers M. In (a) the average values 
(F q ) for the background are all less than, but close to, 1. 
In (b) for the case including 1-jet and in (c) for the bg+2- 
jet case, (F q ) are distinctly higher than 1, especially for 



q > 2. These results represent the first indication that 
the FM method is useful, and that the suppression of 
the statistical fluctuation is realized upon event averag- 
ing without explicit subtraction of the background. That 
is, with (Fq) being around 1 for background only, any sig- 
nificant enhancement of (F q ) above 1 would indicate the 
presence of non-statistical signal. 
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FIG. 1: Event-by-event distributions of F q for (a) q — 2 and 
(b) q = 3 for 2000 events with 30 bins. 
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FIG. 2: Event-averaged (F q ) vs bin number M for (a) bg only, 
(b) bg+lj, and (c) bg+2j. Circles, triangles and squares are 
for q — 2,3, 4, respectively. 

In a more realistic situation elliptic flow introduces 
a tp dependence to the background. We have accord- 
ingly modified the statistical background by a factor 
I + 2u2Cos2iy9. We find that Fig. 2 remains essentially 
the same for v-i increasing up to 0.1 , and that the lj and 
2j cases can be readily distinguished from the background 
with flow. The reason is that the bin-to-bin variation of 
the 2i>2 cos 2tp factor is small compared to the variation 
due to a jet in any event, and the FMs are able to dis- 
tinguish the two types of variation event- by-event. 

In the test cases studied above we have taken the back- 
ground to consist of only 60 particles, distributed over an 
interval of I — 3, corresponding to dN/dtp\b g ~ 20. We 
use dN/ dtp to denote the event-averaged azimuthal distri- 
bution, as is conventionally done, although in our simula- 
tion here the same notation is used for the distribution in 
any event, for brevity's sake. In some experimental data 
the lower bound of p|? soc has been set as low as 0.15 
GeV/c [lj, for which the background height is as high as 
dN/dtp\b g k 200. In such cases the signal of jets on the 
away side is at the 1% level; it is therefore in the presence 
of statistical bin-to-bin fluctuations that can be signifi- 
cantly larger, and our method cannot be expected to be 
effective. We have found that for dN/dtp\b g < 50, the 



3 



lj and 2j signals quantified by (F q ) can be distinguished 
from the background. 

A disadvantage in working with F q is that one no 
longer sees visually the peaks in A4> distribution asso- 
ciated with jets. Thus one gets a quantitative measure of 
the peaks at the expense of losing information about the 
locations of the peaks. However, the loss can be reduced 
if we elevate the level of horizontal analysis of the FM 
by being more specific about the spatial regions in tp. To 
that end let us define 
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where R± stands for the region where ip is < 0, and 
M± = M/2 is the number of bins in R±. Thus we have 
fq = {fq + fq)/2 for each event. We further define the 
NFM, as in Eq. J2J>, by 



Ft = ft /(h) 9 , 



(7) 



which implies F q = (F+ + F~)/2 for every event. 

Since event averaging is likely to erase the distinction 
between F+ and F~ if a jet fluctuates between being on 
the + and — sides, it is useful to consider the difference 
moments D q = |F+ — F~\ and the sum S q = F+ + F~, 
To amplify the effect of the fluctuations let us define the 
vertical average of the pth moment of D q by 



(D p q ) (5) = jf 



JV ovt 



(8) 



and similarly of S q . A plot of (D%) (S) versus (S*) (S) 
for the three cases of bg, lj and 2j for various values of S 
can reveal some characteristics of interest. In our simu- 
lation we place the two jets always on the opposite sides 
of ip = in order to contrast the 2j from the lj events. In 
Fig. 3 we show an array of such plots for p, q = 2, 3, 4. In 
all, the background points are clustered together, while 
the lj and 2j points fan out from the origin, mostly in 
straight lines. Clearly, two jets on two sides of p = re- 
duce D q , resulting in lower (D?) compared to the same 
from lj. Similar plots of the experimental data would in- 
dicate first of all whether the (-D?) and (S?) points of the 
background cluster together at the lower-left corner, and 
then secondly whether the points arising from lj and 2j 
are well separated from those from just the background. 

The implication of Fig. 3 is that the bin size S is not an 
essential variable, and that the normalized asymmetry 



A„ 



(D P q)l(S P q ) 



(9) 



can be a more succinct measure. To see that explicitly 
we show in Fig. 4 the values of A qp vs 5 _1 for the lj and 
2j points. For q = 3 and 4 the dependences on p and S 
are not sensitive, and the differentiation between the lj 
and 2j cases can easily be made, even for q = 2. It seems 
that this type of analysis renders more quantitative and 
distinctive results than 3-particle correlation. 
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FIG. 3: (Dp) vs (S%) for M = 20, 30, 40, 50. 
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FIG. 4: A qp vs S' 1 for (a) q = 2, (b) q = 3, and (c) q = 4. 
Thick, medium, thin lines are for p = 2,3, 4; solid (dashed) 
lines for bg+lj (bg+2j). 



We have generalized the value of Ni to be Gaussian 
distributed, and then also the jet multiplicity m. We 
have found no significant effect on the asymmetry mea- 
sure A qp . Since our aim here is only to suggest useful 
measures, and not their detail numbers, we omit the pre- 
sentation of the results of the Gaussian-distributed mul- 
tiplicities, leaving them to be reported elsewhere. 

Restricting our attention now to only the lj case, we 
consider the question of how best to reveal the position 
and shape of a peak on the away side using FM. To that 
end we need to introduce a cut in ip. Define 

/I (vc) = j^- n ^ - x ) • • • (io) 

< 

where 5> is the set of bins in the range \p\<(p c , and M> 

is the number of bins in S> . The overall FM is then 

< 

given by f q — r</< + r>/>, where r> — M>/M. The 
corresponding NFM is 



Fq<(p c )=fq<(p c )/(h) q 



(11) 
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and the overall NFM is F„ 



r<F< 



>F> The rel- 



ative magnitude of F> vs F^, as ip c is varied, provides 
a measure of the shape of the ip distribution of the jet 
location. Let us define 



B„ 



(12) 
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FIG. 5: (a) Event-averaged (p distribution for J\f e vt = 5000 
with M — 30. The background contributions are not shown; 
their levels are [i] 20, [j] 2, [k] 0.2. The widths and locations of 
the peaks are: [i] (0.4, 0.8), [j] (0.3, 0.4), [k] (0.3, 0). (b) The 
corresponding B4,(<p c ) distributions, (c) dN/dip for [j] and [k] 
with bg included, (d) 1 — B4, for the distributions in (c). 

To test its usefulness we consider three cases for which 
the event-averaged distributions above the background 
are shown in Fig. 5(a); the heights of the background, 
dN/dip\ hg , (not shown) are [i] 20, [j] 2, and [k] 0.2. The 
signal in case [i] is at the 1% level, and corresponds to 
very low p^ soc , whereas case [k] has a sharp peak above a 
low background that can arise from high p^ lg and p^ soc . 
In Fig. 5(b) is shown Bi((p c ) for the three cases, which 
are well separated. In all cases background alone gives B q 
very nearly because either (i^) = (F q ) for ip c > 0.1 
or their sum is much larger than their difference. Note 
that even in case [i] where the double peaks are small 
compared to background, we obtain without background 



subtraction nonzero Bi{ip c ), albeit small. In case [j] with 
dN/ip\i, g = 2 the average bin multiplicity n of the back- 
ground is only 0.2, which is much less than q > 1. Nev- 

> 

ertheless, fq defined in Eq. (10) does not vanish due 

> 

to the fluctuations of the bin multiplicities, so Fq ex- 
ists, though very large. The event-averaged l^F^\ are 

between and 4, resulting in B± to be around 0.5. In 
case [k] the strong peak at ip = results in being 
large compared to F> except when ip c is small; thus B4 
is large, specially at large <p c . Where it deviates from 1 is 
a measure of the narrowness of the peak. B q for q = 2, 3 
give similar results. 

To achieve a better perspective of the mapping be- 
tween dN/dtp and B4, we show in Fig. 5(c) the former 
for cases [j] and [k] with background included. The cor- 
responding curves for 1 — Bi(ip c ) are shown in (d), in 
which [j] exhibits the features of both the broad bump 
in (a) and the high background in (c), while the peak of 
[k] reflects the same in (a) and (c). These properties are 
remarkable due to the drastic difference in the nature of 
the measures displayed. This demonstrates that without 
background subtraction the FM analysis event-by-event 
results in clear and quantitative description of the A<j) 
characteristics. 

In summary, we have considered the use of facto- 
rial moments to analyze the A<j) distribution of parti- 
cles produced opposite a trigger. The advantage of such 
a method of analysis is that no explicit background sub- 
traction is necessary. The FM of each event is sensitive to 
the jet characteristics, while suppressing the effect of sta- 
tistical fluctuation of the background; it is also insensitive 
to the smooth variation due to elliptic flow. The asym- 
metry moments A qp can well separate one-jet and two-jet 
recoil scenarios, and B q can give a quantitative descrip- 
tion of the single-jet characteristics. Application of this 
method to the analysis of the RHIC data on jet corre- 
lation may provide a common framework to compare re- 
sults from widely different experimental conditions and 
various subtraction schemes. 
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